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1 .  INTRODUCTION 


Andreas  (1988)  in  referencing  Wyngaard  et  al.  (1971)  states  that  the  structure 
parameter  c\  should  have  a  universal  form  when  properly  scaled,  that  is, 


cl  z 


2/3 


tt 


(1) 


\ 

Where,  for  unstable  conditions. 


(2) 


z  -  height 

t*  -  temperature  scaling  length 
L  -  Obkuhov  (1946)  length 

Andreas  also  notes  that  Fairall  et  al .  (1980a)  and  Kohsiek  (1982)  extended  these 
arguments  to  the  humidity  structure  parameter  Cp  and  to  the  temperature-humidity 
structure  parameter  C^q  obtaining 


and 


cl  z2/3 


(3) 


(4) 


where , 


Q  -  absolute  humidity  (grams  per  cubic  meter) 

Q*  -  absolute  humidity  scaling  length 

Andreas  further  suggests  that  should  exhibit  the  characteristics  of  ,  Cp, 
and  since  it  can  be  shown  that 

cl  ^  A^CI  *2AB  C^Q  *  cl  <5) 
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where  A  and  B  are  derived  from  a  real  index  of  refraction  fluctuation  expression 
that  varies  with  wavelength.  Consistent  with  the  above  formulations  and 
concepts,  Andreas  developed  a  similarity  expression  for  the  real  index  of 
refraction  fluctuations, 


n.  '  A  t.  +  S  Q, 


(6) 


and  concludes  that 


should  have  the  form 


z2/3 

n! 


(7) 


Later  in  his  report,  Andreas  (1988)  concludes  and  Hill  (1989)  proves  that 

[t]  ‘  [l]  ‘  (t)  *  ^  [t]  *  ^  [t] 


(8) 


We  wish  to  emphasize  that  up  to  this  point  in  this  report  that  Q  represents 
moisture  in  absolute  humidity  form;  whereas,  in  later  developments  we  will  use 
q  (specific  humidity)  to  be  constant  with  Tatarski's  (1961)  original  derivations 

of  C^. 

The  purpose  of  this  report  is  to  present  our  development  of  expressions  for 

('t  >  ^tq  >  >  St  [t]’  [t]’  .  and  (x]  ;  comparisons  with  those  in 

Andreas;  and  comparisons  of  with  Wesely  (1976). 

2.  RACHELE/TUNICK  (RT)  FORMULATIONS 

We  initiate  our  development  using  Tatarski's  expression  for  ,  that  is. 


1*1  fer 


(9) 


where , 

b  -  constant  -  3.2  (Wyngaard  (1973),  Hill  (1989),  Andreas  (1988)) 

Kg  =  -  ^  *  turbulent  exchange  coefficient  for  heat  (10) 

9h 
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energy  dissipation  rate 


(11) 


I  <  0 


(12) 


L  <  0 


(13) 


k  —  Karman's  constant  (0.4) 
u*  -  friction  velocity 
7  =  ^  =  15  (Hansen.  1980) 

^  “  the  gradient  of  the  real  index  of  refraction  in  the  vertical 
dz 

direction  z. 

Andreas'  (1988)  expression  for  n  with  moisture  in  absolute  humidity  is 


n  ^  A  (X,  P,  T,  Q)  T  *  B  (A,  P,  T,  Q)  Q, 


(14) 


where  for  X  =  0.55  pm 

A  -  -  7.9  X  10-5  JP 
B  -  -  5.46  X  10-5 

and  where  P  and  T  are  pressure  in  millibars  and  temperature  in  K. 
In  our  notation  we  write  equation  (14)  as 

dn  ‘  A  (X,  P,  T,  q)  dT  *  B,  (A,  P.  T,  q)  dq 


(15) 


where, 
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A  -  -  7.9  X  10-5 


B,  =  -  1.97  X  10 


-5  P 


q  -  specific  humidity 
Dividing  equation  (15)  by  dz  gives 


^  =  -  7.9  X  10-5  P  ^  -  1.97  X  10-5  P  ^ 
^  72  ^  T  dz 


Furthermore,  equation  (5),  in  terms  of  specific  humidity,  becomes 


CI^IAB  C^.  +  cl. 


For  unstable  homogeneous  conditions,  we  let 


i  and 

Si  Si  Si 


Substituting  equation  (18)  into  equation  (16)  gives 


=  -  7.9  X  10-5_L  4^  +  -  1.97  X  10-5  4  4^ 
Si  Sz  T  Si 


^  ~  -  7.9  X  10-5  P  4«  -  (1.975  X  10-5r)  4^ 

Si  Y^Sz  Y^  Si 


Squaring  equation  (20)  gives 
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dn 

■3i 

.  j 


=  6.241  X  10-® 


p2 

2»4 


+  [3 . 11  X  10-®  r]  ll 


de  aq 
Tz  ~3z 


+  [3.88  X  10-1°  1.2] 


p2 

Jl 


dO  do 

Expressing  and  in  equation  (21)  In  similarity  form,  for  unstable 
ditions,  that  is, 


88 

■3i 


equation  (21)  becomes 


dq 

■31 


=  6.241  X  10-«  £! 

r'*  Jc^ 


..2 


+ 


3.11  X  10-8  j 


p2 


+  3.88  X  10-1°  j-2 


El 

r* 


2 


q 


where 


7  -  constant 

L  =  U*'  Til  ■*-  0.61  Q) 

kg  (0*  +  0 .61  0  g* 


(21) 


con- 


(22) 


(23) 


(24) 


(25) 
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g  -  acceleration  of  gravity. 


Substituting  equation  (2A)  into  equation  (5)  and  (9)  gives 

Cn  =  A'  +  Z'  e'  q*  +  C'  q*^ 


(26) 


where 


6 . 241  X  10-®  ^  .  (l 


(3.88  X  10*1°  r2) 


(27) 


(28) 


"?9) 


€  = 


k  u*  z 


(30) 


(31) 


(32) 


From  equations  (30)  through  (32) 


(33) 


Substituting  equations  (33)  into  equations  (27),  (28),  and  (29)  gives 
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(34) 


A'  =  b  (6.241  X  10-3  Ij-2/3  ^-2/3 

=  b  [3.11  X  10-3  j]  ._Pf  ^-2/3  ^-2/3 

C'  =  b  (3.88  X  10-^°  r2)  ^  k-2/3  z-2/3  ^1  -  7  II 
Equating  equations  (26)  and  (17)  term  by  term  gives 


(35) 


(36) 


(37) 


(38) 


(39) 


(40) 


(41) 
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b  rq-ic-2/3|  |  = 


or 


C  z2/3 

tq  ^  ^  h  t-2/3,  V  ,  „ 

e~q-  W  «tq 


{}■ 


Also 


=  b 


iC-2/3 


or 


C,"  ^2/3 

"7^ 


6  ic'2/3 


To  account  for  moisture  contributions  to  ,  Kunkel  and  Walters  (1983) 
expression  by  Wesely  (1976),  that  is, 

Cn  =  cl  a2  (1  +  0.03/Bo)2 


where 


Bg  =  Bowen  ratio 


4  X  10-*il 

q* 


In  our  model  we  relate  to  and  as  follows.  From  equations  (37) 
(45)  we  write 


1.2 


2-2/3 


8t 


.42) 

(43) 

(44) 

(45) 

used  an 

(46) 

(47) 

through 

(48) 
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'tq 


S,,  =  r  q-  z-"/'  gt 

(49) 

gc  = 

(50) 

Furthermore,  from  equations  (48)  through  (50),  we  obtain 


r  -  <1*  Stq  _  q*  ^2 

«  ■  r  ‘  X  7"  ‘ 


^2  _  q  Sq  y^2 


(^r 


Substituting  equations  (51)  and  (52)  into  equation  (17)  gives 

ci 


31*  isT  tiiT 

A  r  l^J  [rJJ 


(51) 


(52) 


(53) 


or,  in  terms  of  Bq, 


T  ^  Bq  ^  1.6  X  10-^  fBqf 

ab;  — 


(54) 


Substituting  the  values  of  A  and  B^,  given  earlier  (equation  (15)),  into  equation 
(54)  results  in 


Cl  [l  .  1-99x10-*  j.  ,  3.98x  l0-«  j2 


B‘ 


(55) 


3 .  COMPARISONS 

Wyngaard  (1973)  concluded  from  scaling  laws  and  the  analyses  of  Kansas  data  that 

g,(0 ' 4.9  (1  -  6.1r)-2/3  (56) 
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for  r  ^  0.  Hill  (1989)  demonstrated  that  '  gtq(r)  =  gqCD  =  • 

Andreas  (1988)  concluded  that 


gb(r)  = 


5.92^H(f) 

■TTcrF^ 


(57) 


based  on  the  budget  equations  where  Wyngaard  and  Cote  (1971),  Kaimal  et  al . 
(1972),  Fairall  et  al .  (1980b),  and  Schacher  et  al .  (1981)  found  that 


4>e(n  =  [l  ♦  0.46(~O2/3f'^  2  <  f  <  0 


^.(O  =  [l  +  2.3  £-3/5  0  <  C  <  2 


(59) 


From  the  formulations  in  this  study  (equation  (40)),  we  found  that 


gt(r)  = 


5.89^g(n 


(60) 


where,  from  Panofsky  (1968),  *=  ((|>„  ~ 

We  also  show  that  gt,(0  =  gq(0  =gtq(C)  =Sn(^)- 

The  functions  (1  -  6.10'^^^. _ ^j-g  compared  in 

[l  >  0.46(- 

figure  1;  whereas,  the  complete  functions 

4.9(1  -  6.  If)  "3/3^ - L'  ^ - _,  and  ^  f  are  shown  in  figure  2. 

[l  *  0.46(-  '  (^n.  - 

Figure  3  illustrates  resultant  vertical  profiles  by  using  the  three  forms  for 
g(z/L)  discussed  above.  Equations  (6)  and  (7)  describe  the  formulation  for  the 
profiles  presented  in  figure  3. 

Finally,  we  compare  values  of  C^/Cl  by  using  equations  (46)  and  (55).  These 
values  are  presented  in  figure  4. 
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V' 

^  V' 

^  V 

\NX 

Rochele. - 3oid  ^ 

Andheos. - lg.dashed  ^ 

\Aytgoard — .smdaahed 


lor*  W’ 

-2/L 


Figure  1.  A  comparison  of  the  variations  of  the 

similarity  portion  of  the  function  g  (z/L) 


Figure  2.  A  comparison  of  the  complete  expressions 
of  the  function  g  (z/L) . 
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S-|  RcxMe. - aotd 

Andeos. — Jg.dashed 

^  wyn90<rt.-_*xd«»»d  5^ 

f~~  I  I  I  I  1 1  m  I  I  1  n  1'  I  M  f 

XT*  10"**  ior“ 

Ch^ 

Figure  3.  vertical  profiles  using  the  g  (z/L)  functions  of 

Rachele,  Andreas,  and  Wyngaard. 


Weseiy. - dashed 


Figure  4.  The  ratio  C^/C^  versus  Bp  (Bowen  Ratio)  for  equations  (46) 
(Wesely,  1976)  and  (55)  (Rachele)  for  T  -  296.15  K. 


4.  SUMMARY  AND  CONCLUSIONS 


In  this  report,  we  have  shown  our  development  of  expressions  for  the  refractive 

index  structure  parameters,  ,  as  well  as  the  temperature  and  humidity 

structure  function  parameters  ,  and  C^g.  We  have  also  presented  our 

development  for  the  similarity  functions  gj^,  and  gg. 

We  concluded  that  the  singular  universal  function  g  (z/L)  satisfies  the 
expression  gt  -  gtq  -  gg  -  g  as  shown  by  equations  (40),  (43),  and  (45). 

In  figures  1  and  2  we  show  that  as  I  "-j  I  increases  in  magnitude  greater  dif¬ 

ferences  arose  between  the  expressions  for  g^  of  Andreas  (1988),  Wyngaard  (1973), 
and  Rachele.  Clearly  they  are  not  equivalent. 

However,  much  less  separation  is  evidenced  between  the  profiles  of  pre¬ 
sented  in  figure  3.  The  choice  of  which  g(z/L)  function  to  use  may  now  be  less 

critical . 

Finally,  we  compared  values  of  the  ratio  C^/C^  by  using  equations  from  our 
development  and  those  from  Wesely  (1976).  Figure  4  shows  that  both  expressions 
compare  closely  to  one  another. 
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